Nonlinear Dynamics of Quasiclassical Spin Moment in a Swept Field 
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Quantum dynamics of anisotropic spin system with large spin moment in a swept magnetic field 
is theoretically investigated. Magnetic field of this type induces vortex static electrical field, that 
breaks down the axial symmetry and induces new coherent quantum phenomena: the occurrence 
of band energy spectrum with continuous spin states, Bloch-type oscillations and interband Zener 
tunnelling effect. These quantum phenomena display themselves, in particular, in magnetization 
jumps and susceptibility peaks in the investigated spin system. 
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1. The problems concerning to quantum dynamics of 
spin systems have aroused considerable interest in recent 
years (see references § ||, |, | §). In many respects it 
is interlinked with recent discoveries of the macroscopic 
quantum tunnelling of magnetization Jl], ^[ ||, [I| ^] , quan- 
tum hysteresis, molecular bystability, novel type of mag- 
netic oscillations associated with Berri phase appearance. 
These mesoscopic phenomena were discovered in so called 
"giant spin moment systems" - in the magnetic nanoclus- 
ters Mni2, Fes with large spin moment in the ground 
state (S = 10). Some mesoscopic spin phenomena are 
likely to be discovered in the rare-earth ions with large 
angular momentum (Dy 3+ , Tb 3+ , Ho 3+ , Er 3+ ). Of par- 
ticular interest are problems pertaining to the macro- 
scopic quantum coherence, quantum measurements in 
spin systems and the quantum correlations violation due 
to interaction with environment and especially quantum 
to classical behavior transition. It should be stressed 
that the observing conditions of these coherent phenom- 
ena are more strict than that for above-listed. Therefore 
it is seemed to be of importance to explore the novel 
situations when quantum coherence phenomena appear. 
This problem is of practical concern for magnetic nano- 
electronics (spintronics) and quantum computer science. 
One of the most intriguing ideas in spintronics is the 
"giant-spin" -nanocluster usage as a bistable unit for fur- 
ther generation of the molecular memory. These nan- 
oclusters may be interesting as a promising realization of 
qubits in quantum computing ]?], ||, ||, [nj . 

The aim of the paper is quantum dynamics investi- 
gation of the large spin moment quantum system in a 
swept magnetic field. Such a field produces a torque on 
the spin system inducing the spin precession and thus, 
displays new effects in the dynamics of the spin system. 
The present paper develops the ideas surveyed by author 
earlier as applied to the antiferromagnetic nanoclusters 
Jj"l| , metal rings and ring-like molecules [[l2| . 

2. We consider a quantum system (ion, molecule, clus- 
ter) under the swept magnetic field. The Hamiltonian of 
the system is taken to be 



H = gns JB(t)+V CF , 



(1) 



namics of the investigated system the quasiclassical ap- 
proach will be used. 

The crystalline field is represented as 

Vcf = V CF + V CF , 

where Vq F is the crystalline field with "easy plane" sym- 
metry and V CF creates the anisotropy in plane; in so 
doing \ Vl; F \ <C \ Vq F \. Let z axis of the Cartesian coordi- 
nate system be perpendicular to the "easy plane" which 
contains 2-nd order axis selected as x-axis. 1 

Assume that the magnetic field is B = (0, 0, B z ) with 
B z being time-dependent. Let Bit) = Bit/r, where B\ 
and t are the characteristics of the field increase (de- 
crease) process. According to the well-known analogy let 
designate the j m = B\ /4ttt as a "magnetic current" . 

3. To describe the spin dynamics the coherent quan- 
tum states technique \9, ip) Jjl| [l4| will be employed with 
6 and if being polar and azimuthal angles of the angular 
momentum. These angles are measured from the z and 
x axes respectively. 

The Lagrangian of the system can be represented as: 



M 

C = {cos9 - 1) <p - E, 

1 



E = -K 1 sin 2 6 + K 2 sin 2 9 sin 2 p - MB(t) cost 
|tf 2 | « #i>0, B\\z. 



(2) 



(3) 



This Lagrangian can be deduced with conventional co- 
herent quantum states technique iQ. All terms in (Q) 
are of apparent physical meaning. The first and the sec- 
ond terms in (|) are: the so called Wess-Zumino term 
(kinetic energy) appeared because of the nonorthogonal- 
ity of the coherent states at t and t + At and the total 
energy E. The first and the second terms in (^) are the 
anisotropy energy (crystalline field). The 8 and tp func- 
tional dependences of E are chosen in the simplest but 
sufficient way to elucidate the principle aspects of the 



where Vcf is the crystalline field operator. It is suggested 
that J > 1, thus in order to describe the quantum dy- 



The generalization for situations with the anisotropy in plane 
being characterized by other symmetry elements, for example 
4 Z , 3 Z , 6z, is not a particular problem. 
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problem. The last term in (Q) is Zeeman energy. The 
Euler-Lagrange equations for Lagrangian (||) are equiva- 
lent to the Landau-Liphshits equations (without attenu- 
ation). 

4. The partition function of the quantum system can 
be represented as the functional integral in the Eucledian 
space (t = it) 



Z= DcosO / D<peif° 



(4) 



where = ±, = 9(t), p = ip{r). As \K 2 \ < |JCi| 
one can suppose that in moderately high magnetic fields 
9 — 7r/2 <C 1. Then integrating the partition function 
expression (Q) upon 9 yields 2 : 



Dip e h J o 



drC e 



(5) 



where 



£e// = -I£--U A (<p)->yI + iB(t)<p, (6) 



where / = m r 

U A = Ki sin 2 <p. 



is the moment of inertia of the system, 
The last term in (H) is defined as the 



Zeeman energy contribution. 

The Lagrangian (^|) can be presented (in the real time) 
as following: 



-eff 



I(<p->yBf 



K2 sin ip 



(7) 



It should be specially noted that the variable ip is de- 
fined here not on the S 1 (0 ^ p < 2ir) set as it usu- 
ally assumed in the theory of angular momentum but 
on the real numbers V set (p € R 1 ). The latter in the 
present problem represents a product bundle of the S 1 
set. S* 1 plays a role of the base of this fiber set V. This 
remark is obviously important in the given context as 
presence of the swept field B z (t) breaks down the sym- 
metry ip — * p + 2im transformation, where n is integer, 
therefore, S 1 set, commonly used in such problems, must 
be enlarged up to the V ■ The axial symmetry violation is 
apparent because the swept field B z (t) generates the elec- 
trical axially symmetrical field E v , by virtue of Maxwell 

equation — ~B — curlE. 

Systems with potential energy of the "tilted washing 
board " type U(x) = Uq(x) + cx, where Uq(x) is the pe- 
riodical function of x, c is constant value, were already 
investigated earlier. As for example it is possible to men- 
tion the electron motion in electrical field in a crystal 
|pp| , [2l|| or Josephson junction dynamics when direct cur- 
rent flowing (TJ, |T7|. Therefore one can expect the man- 
ifestation of some properties in the spin momentum dy- 
namics which are similar to that of above discussed sys- 
tems. Such characteristic properties are the band energy 
spectrum occurrence, the Bloch oscillatio ns |20|, pl[ and 
the interband Zener tunnelling effect |^2|, [23^ |24|. Let's 
consider these issues in more detail. 

5. The generalized momentum appropriate to the co- 
ordinate ip equals to p v = q^ s — J ip — I(ip — jB). 
Then the Hamiltonian of the system TL — p v p — C e ff 
can be represented as 



which is convenient for the gauge transformation per- 
formance and for making analogies with other quantum 
systems (the Josephson junction |l5|, fl6[ , m esoscopic con- 
ducting ring or ring- like molecule |12|, |17| , antiferromag- 
netic nanocluster One can readily see that the La- 

grangian ([?]) can be recast with an accuracy of the total 
time derivative as (see also 19 ) 



C 



eff 



2 



Kn 



cos 2p + jIBip. 



(8) 



where P v — 



H=±(P v - 7 IBf + U A (<p), 
h d 

i dip ' 



(9) 



The gauge transformation ^i(ip) — ► ^(f) e 9 ^'^, where 
9 = j:I~/B p (the label "1" of the function * will be 
further omitted), transfers the Schrocdinger equation to 
the following type: 



Pi 
21 



+ U a (cp)-IjB<p 



(10) 



2 The calculation procedure can be reduced to the following. It is 
not difficult to see that value 8 = ir/2 is the extremum point of 
the Eucledian action of the functional Se = Jq^ dr C(9,tp) (at 
B = 0). From the stationarity action condition (Saddle-point) 
SS E = it follows (at B + 0): 



■0 M 



B(t) 



(4') 



M 

2 "' ~2Kl \ " 7 

Substituting (Ej) to (^), then expanding S E about A6 = (8-8 M ) 
to second order, and calculating, if need be, the Gaussian integral 
upon A9, equations feh and (6) can be obtained. In equation (|5h 
the Gaussian integrafis omitted for simplification of the formula 
because of being further unessential. 



At first let's consider the eigenvalue problem of the 
Hamiltonian (|j|) at B = . The eigenstates of the Hamil- 
tonian (^) are the Bloch functions 



y s (p + 2n) = e™ m M> s (p), 



(11) 



where m is an arbitrary real number, s is the energy 
band number. It would appear reasonable to designate 
the parameter m as quasispin, (compare with the quasi- 
momentum for a band electron). By analogy with the 
term - "charge states"- for similar states in the theory 
of Josephson effect it is possible to call the states (|T^) as 
the "continuous spin states" . It is common knowledge for 
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the spin moment component to be quantized at the fixed 
direction. In the case being considered the "quasispin" 
is a arbitrary real number (m G R 1 ) 

The difference between these two types of states can be 
explained as follows. The quantized spin states are spec- 
ified at 5* 1 set (0 ^ ip < 2ir) , in so doing the spin moment 
quantization is clearly associated with the axial symme- 
try of the system, in other words, with the boundary 
conditions , 5(<yS + 27r) = ±^((p). The absence of this sym- 
metry in the dynamic symmetry group of the Lagrangian 
(||) cancels the quantization of a spin moment. Instead 
of this the "continuous spin states" i.e. Bloch functions 
(|TT|), are achieved in the fiber set V (oo < ip < oo). 

Let Ua(<p) — — (1/2) 0082(^9, where K 2 is a con- 
stant. Then the Schroedinger equation for the Hamil- 
tonian (||) is reduced to the Mathieu equation from the 
theory of which it follows that the energy spectrum of the 
Hamiltonian (0) has a band structure, i.e. the eigenvalues 
of (^|) E n (m) are the functions defined in the appropriate 
Brillouin zones. At K 2 ~ the band structure corre- 
sponds to the approximation of "free electrons" type: 



E s (m) 



h 2 m 2 
21 



(12) 



with the forbidden bands at the boundaries of the Bril- 
louin zones: mg = s(s = ... — 2, —1, 0, 1, 2, ... ) which 
are narrow in accordance with K 2 <C K\ . 

Near the Brillouin zone boundary, for example, near 
m b = —1, the wave function can be represented as fol- 
lows 



-iEt 



(13) 



where u(ip) = A ie mip + A 2 e l{m+2)ip . The Schroedinger 
equation for the Hamiltonian (^) can be written as (the 
Mathieu equation) 



u" + (fi — 2d cos 20) u = 0, 



where [i 



2IE 



d = 



IKo 



(14) 



2h 2 ' 

<p + 7r/2. The sign "<- 



Used here is a new variable <p 
will be further omitted. Substituting ( |l3| ) to (|l4| ) yields 

(/i - m 2 ) Ai - dA 2 = 0, 

- dAi + (/i - (m + 2) 2 ) A 2 = 0, 

whence it follows that 



(15) 



m 2 + (m + 2) 2 
A« = I ^ — ± 



(m 2 - (m + 2) 2 Y 



d 2 



(16) 



In particular from the expression (|16|) it is clear that at 
m.B = — 1 the forbidden band width (i.e. between allowed 
zones with s = and s = 1) equals 



AE 



2h 2 d 



10 



= K 2 



(17) 



The allowed zone widths of the zero, first and third zones 
are E « E x w ^ and E 2 w ^f-, respectively, 
whereas the forbidden band widths rapidly decrease (at 
K 2 <C K\) when going to the higher energy bands. Then 



A£ 2 i 



AE; 



32 



AE d 



16/ ^2^ 2 
ft 2 /7^ 2 



IK2 

2h 2 

IKo 



(18) 



576/ V 2^ 2 



etc. 



The equations @@, @, @ determine with suffi- 
cient accuracy the energy spectrum in the limit of first 
two Brillouin zones. 

6. Let's discuss effects of the swept magnetic field 
B z (t) effects that as mentioned above can be considered 
as a classical gauge field in equations (|9|,|l0|). The last 
term jIBip in equation (^) is playing the same role as 
the energy - eF-x (F characterizes the electrical field and 
x is the electron coordinate) in the well-known problem of 
the dynamics of the Bloch electron in an electrical field. 

Let consider the momentum p v dynamics in the case 
of magnetic field varying adiabatically slowly: 



J7B «Jf 2 . 



(19) 



To describe the dynamics of a spin under the action of 
"magnetic current" j m = B/Att let's consider a wave 
packet of the Bloch functions ([ill). Let m and (p be the 
mean values of the quasispin and the coordinate of the 
wave packet center and the values Am, Aip (Am-Atp ~ 1) 
determine the corresponding uncertainties. Under influ- 
ence of "magnetic current" j m the wave packet formed 
at t = moves to the boundary (for example, right, i.e. 
nig = 1) of the Brillouin zone where Bragg reflection 
takes place. Here the wave packet velocity changes its 
sign to an opposite and the wave packet appears at the 
left boundary {tub = — 1), after that it moves again to 
the right boundary of Brillouin zone (mg = 1) then re- 
flects again and so on. This process is called the Bloch 
oscillations. It is described by the following equations for 
m and (p mean values: 



Tit 



1 dE s {rh) 



(20) 



dm 



During this (adiabatic) process the system remains in the 
state with the definite s and its physical properties as, 
for example, the magnetic moment, are the oscillating 
functions of time with the frequency 



Bloch — 



(21) 
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If the external magnetic field has the harmonic contribu- 
tion i.e. 



B = Bit/r + w sin27r ft, 



(22) 



then the resonances on the frequences / = f Bloch and 
/ = r f Block, where r is a rational number, are available 
(the Stark ladder-like resonances). 
When the "magnetic current" increases 



jIB 



>K, 



(23) 



the Zener tunneling effect occurs between adjacent zones. 
In particular, the probability of the tunnel transition in 
a unit time between zones with s — and s = 1 equals 



.9oi — f Block e 



(24) 



where 8 — Zf 2 J , 7 = — . 

" H *yBl 5 ' rac 

7. Let us consider the average magnetic moment be- 
haviour of the considered spin system qualitatively. The 
magnetic moment z-axis component is 



M 7 = M cos 9 « M 



(§ 



M 



f2 



2Ki 



B, 



(25) 



Averaging (]2 
yields 



with the corresponding wave function 



(M, 



B, 



1 



(26) 



where x± = jk^- 

Two limiting cases have to be considered: a) K2 = 
therewith g n ,n±i — 1 and b) goi = 0. The first case can 
be defined as a free precession of the spin system and the 
second one is a Bloch oscillations case. 
At K2 — expressions (Eoh 



give 



(<p)=j(B(t)+c) 



(27) 



where c is a constant value determined by initial condi- 
tions. Substituting @ to @ yields 



(28) 



(M) = M , 



where Mq = — \±_ c is the system magnetic moment at 
t = 0. Thus, the distinctive property of the free preces- 
sion case is that the accelerated spin precession influenced 
by increasing (decreasing) magnetic field screens the ion 
paramagnetic susceptibility contribution (xj_B z ) so the 
average magnetic moment is field undependent. 

In the case of the Bloch oscillations the situation 
changes drastically. Here the (M z ) versus B z dependence 
represents the sum of "ordinary" linear (x± B z ) term and 
the periodical curve (ip) with the period of 



AB = Bi (t f Block) =%hl. 



(29) 



At Ki <C K\ this periodical function is rather close to 
"Saw- type" with the amplitude Aip ss 21 /H: 



B, sC B < b, 

B -2b, b^ B <3b, 

B — 46, 3b < B < 5b 
etc., 



(30) 



where b = \. 

fl 

Generally, the M Z (B Z ) dependence contains the char- 
acteristic features of both limiting processes. 
Actually, as it was shown above (see (|l6|)) at K-xjK\ <c 1 
the forbidden band AE^i width falls far short of AE±o, 
AE32 *C AE21 etc. Therefore as a first approximation 
it is possible to put 512 = .923 = <?34 = •■• = 1, i-C- 
to neglect the Bloch oscillations in the first and follow- 
ing excited bands. It means that the precession in all 
bands except ground one (i.e. s = 0) can be considered 
as free. The precession influenced by B z (t) as a whole 
can be assumed as following: the wave packet formed at 
t = slightly spreading out reaches the edge of a Bril- 
louin zone, then partially reflects with probability 1 —.901 
and partially tunnels (with probability goi) to the next 
zone where it freely precessing. With this considerations 
and the equations (p4|),(|30|) the M Z (B) dependence can 
be presented as follows: 



0, 

26(1 -p), 
Xx X M z = <(4&(l-p) 2 , 
66(1 -pf, 
etc., 

V 

- „p 



< B < b, 
b < B < 36, 
36 ^ B 56, 
56 ^ B € lb 



(31) 



where p — e" , (3 is determined in (g4j). It is easy to 
verify that the magnetic moment step AM Z magnitudes 
decrease as the number of the step increases. So, at B = b 
(AM Z ) W = x± 26(l-p); at B = 36 (AM,) 2 i = X± 26(1- 
p)(l - 2p); at B = 56 (AAf z ) 32 = x± 26(1 - p) 2 (l - 3p) 
etc. This means that the M Z (B) dependence tends to 
saturation when i? increases. 

All above mentioned considerations are concerning to 
the case of T = K. Obviously, the thermal fluctuations 
(at T ^ OK) and dissipation (i.e. dissipative environment 
interaction) breaks down the studied quantum coherent 
effects. Taking into account the case of finite tempera- 
tures is worth another look. Here we shall restrict by only 
pointing out that the "T = K" -theory is valid when 
T <C K\ and T 2 3> br/Bi, where T2 is spin relaxation 
time of the system. At K\ ssl cm~ 13 and B\/t S3 10 4 



3 It is worth to note that equation looks realistic for a triplet 
ground state originating from the crystalline field splitting of the 
ground multiple of the spin system. 
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Oe/s we have T <2 K, T 2 »0.04 s. These limitations 
appears as quite executable for modern low temperature 
experiments. 

It should be stressed that the neglected above Bloch 
oscillations in the bands s = 1, s = 2, etc. and the Bragg 
reflections at the corresponding Brillouin boundaries give 
the additional periodical terms into the M Z (B) and x(-B) 
dependences. So, the Bragg reflections at the boundary 
(s = f, 8 = 2) give the additional intermediate peaks 
of x at the fields B — 26, 46 etc. The reflection at the 
boundary (s = — 1, s = 0) decreases the peak amplitude 
at the field B = 36. The Zener tunnelling at the bound- 
ary (s = — 1, s — 0) even more complicates this picture 
due to an interference of the tunnelling wave function 
(s = — 1 — > s = 0) with the main wave function (s = 0). 
These questions are worth to be further investigated in 
more detail. 

In conclusion, it was shown that the external time de- 
pendent magnetic field induces new coherent quantum ef- 
fects of the anisotropic spin system dynamics. They are: 
the occurrence of band energy spectrum with continuous 
spin states, quasibloch oscillations and interband Zener 
tunnelling effect. These quantum effects manifest them- 
selves as magnetizaion jumps and susceptibility peaks in 
the investigated spin system. 
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